The Instantaneous Bethe-Salpeter Equation and Its Analog: the Breit-like Equation 
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We take (/i ± e =F ) systems and consider the states with quantum number J p = 0~ as examples, 
to explore the different contents of the instantaneous Bethe-Salpeter (BS) equation and its analog, 
Breit equation, by solving them exactly. The results show that the two equations are not equivalent, 
although they are analogous. Furthermore, we point out that the Breit equation contains extra un- 
physical solutions, so it should be abandoned if one wishes to have an accurate description of the 
bound states for the instantaneous interacting binding systems. 



^1- 

o 
o 

•*-> 
O 

O 



> 

o 

&\ 
O 

O 



O 
=5 



X 



PACS numbers: ll.10.St, 36.10.Dr, 12.20.Ds 

Keywords: instantaneous BS equation, Breit equation, exact solutions 



In the reference , we pointed out that the authors of 
Refs.@,El etc replaced the Bethe-Salpeter (BS) equation 
which has an instantaneous kernel with an analog of it: 
Breit equation 5, 6] un- properly. To explore the different 
contents i.e. the un-equivalence of the two equations, 
in the paper we restrict ourselves to the bound states 
of the systems (/i ± e =F ) systems with quantum number 
J p = CP (S-wave) to solve the two equations exactly 
and to examine the obtained solutions accordingly. 

BS equation for a fermion-antifermion system has the 
general formulation |^: 



mi)xp{q){^2 + m 2 ) = i 



d A k 



V(P,q,k) XP (k) , 
(1) 

where xp(q) is the BS wave function, P is the total mo- 
mentum, q is relative momentum, and V(P, q, k) is the 
kernel of the equation, pi,p 2 are the momenta of the con- 
stituents 1 and 2 respectively. The total momentum P 
and the relative momentum q are related to the momenta 
Pi,P2 as follows: pi = ol\P + q, ol\ 
a 2 P 



mi 



mitmj ; 



P2 



mi+m 2 ' 



If the kernel V(P, k, q) of the BS equation has the be- 
havior: V(P,q,k)\p_ Q = V(q,k) (P = in center mass 
frame of the concerned bound state) , the BS equation is 
called as an 'instantaneous BS equation'. 

For instance, in Coulomb gauge, the terms correspond- 
ing to the possible transverse-photon exchange between 
the two components in BS kernel for (/i ± e =F ) systems are 
considered as higher order, so the 'lowest order' BS equa- 
tion kernel for the systems has the form: 



V(P,q,k)\ post 
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Namely, there are some physical systems whose bound 
states are described by an instantaneous BS equation in- 
deed. 
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In Refs.pi 0, we showed how to solve an instanta- 
neous BS equation exactly, and also showed the authors 
of Ref.0, HI how to mislead the problem to its analog: 
Breit equation. For present convenience, let us repeat 
the main procedure in Refs.0,0]- 

Firstly, as done by the authors of Refs.0, Q the 'in- 
stantaneous BS wave function' <pp(q) as 



ipp(q) = i J ^-xp(q°,q), 



(3) 



is introduced, then the BS equation Eq.QJ can be re- 
written as 



Xp(q l> ,q)=Sf\p^r ] (q)Sy'(-p^ 



(2), 



(4) 



Here S^\pi) and S?\— p 2 ) are the propagators of the 
fermion and anti-fermion respectively and 

d 3 k 



v(q) = 



(27T) 



V(q,k)vp(k). 



The propagators (in C.M.S. i.e. P = 0) can be decom- 
posed as: 



uo f y.jp l ) - Jq0+a . M _ 



Kit) 



A" (g) 



(5) 



with = y/mf + q 2 , A i (q) = 7% ± J (m 2 + 7 • 

q) , where J = 1 for the fermion (i = 1) and J = — 1 for 
the anti-fermion (i = 2). It is easy to check 

Af(q) + A?(q)= 7 °, 
Af(q) 1 °Af(q) = Q, At(f)~f°At(q) = At(<f). (6) 

Thus A^ can be considered as 'energy' projection oper- 
ators, and 'complete' for the projection. For below dis- 
cussions let us introduce the notations <fp^(q) as: 



^(q) ee A±(g) 7 Vp(9l7°Af (q) 



(7) 



Because of the completeness of the projection for A ± , we 
have: 

(fip(q) = ^ + {q) + (Pp~{q) + f P + {q) + <Pp~{q) 
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for the BS wave function (fp(q). If we further carry out 
a contour integration for the time-component q° on both 
sides of Eq.(@J, then we obtain: 



(M — ZJ\ — UJ2 ) 



(M+uj 



(M is the eigenvalue). To apply the complete set of the 
projection operators Af p (q) to Eq.(jHJ, we then obtain the 
four equations: 

(M - «i - io 2 ) V + P + {q) - A+ (9)7^ , (9) 



{M + lo 1 +lo 2 ) Vp ($) = -Ai(q)r)p(q)Az{q), (10) 



(ii) 



Note that the two equations in Ea. l|lll) do not contain 
eigenvalue M, so essentially they are constraints for the 
BS wave function ipp(q). Because of the completeness of 
the projectors A ± , only the four equations contained in 
Eqs.£3E3inj are equivalent to the equation Eq.©. 

E.E. Salpeter [2| and the authors in literature such as 
Ref. |3 would like to connect the above coupled equations 
to the Breit equation [5j,|6|], so they did a 'combination' 
for Eas. (|9llUlll|l into one operator equation (in C.M.S. 
of the bound state, P = 0): 

[M — Hi(q) — H 2 (q)](p(q) 
= A+(q)^ V (q)^A+(q) - A^ (q) 7 °v(q)l°A^ (q) , 
Hi(q) = 7Ti i/3 + q- a , H 2 (q) = m 2 (3 - q ■ a , (12) 

with the definitions (3 = 7 , a = (3<y. In fact, the 
equation Ea. (|12|) is not equivalent to the instan- 
taneous BS equation i.e. the coupled-equations 
Eqs.lJUl 1101 II 111 . When applying the project oper- 
ator A^(q)7° ® j°A + (q) to Eq.lO we obtain Eq.@, 
when applying the project operator Aj~(g)7° ® 7°A^~(<t) 
to Ea. (|12f) we obtain the Ea. (|lU|l . whereas when apply- 
ing Af (q)j° ® 7°Aj(q) to Ea.lJTIZ)). then we obtain the 
homogeneous equations: 



[M - CJi(q) + UJ 2 (q)]>fi + - (q) = , 
[M + uj 1 {q)-u 2 {q)]^-+{q)=0, 



(13) 



with Wi 2 



q 1 . The two equations in Ea. H13|l 



contain the eigenvalue M, so they do not play a role as 
constraints any longer as played by Ea.l(TT)l. They are 
homogeneous equations so we may conclude that they 
formally not only have the 'trivial solutions' i.e. Ea. l|ll|) . 
but also have 'non-trivial solutions'. Note that when 
mi = m 2 also means u>x = to 2 , so from Ea. (|13|l formally 
one may be sure that only M = will correspond to the 
'non-trivial solutions' of Ea. (|13|) . whereas, when M = 0, 
there is no the frame P = at all, so Ecis. I|9I10I11I12I13|) 
written in P = system will not be correct and we have 
to re-derive them. To avoid the mi = 7712 case, below 



we restrict ourselves to consider the cases mi / 7772 only, 
and it is why we take (/i ± e =F ) systems as examples below. 

It is interesting to explore the differences of the two 
equations: the instantaneous BS equation Ecis. H9I1UI11|) 
and the Breit one Ea. (|12l) or say (|9I10I13I) . because it is 
not only more precisely to prove the un-equivalence of the 
two equations than the above argument about the homo- 
geneous equations Ea. (|13fl . but also is useful for under- 
standing the two equations: the instantaneous BS equa- 
tion and the Breit equation deeply. Specifically, to do it, 
we restrict ourselves on the (/i ± e =F ) systems for the bound 
states with quantum number J p = CP as examples, and 
solve the coupled equations Eas. H9ll0lll|) corresponding 
to the instantaneous BS equation and Eas. (j9ll()ll3|) cor- 
responding to the Breit equation respectively without 
any approximations. In fact, for the present purpose at 
the moment, it is enough only to consider the bound 
states with quantum number J p = CP . 

As pointed out in Ref. pj , in general, the wave functions 
for J p = CP states in C.M.S. (P = 0) have the following 
formulation: 



<t> = (7° '91 +92 +^93+^1° 'gi) 7 5 



(14) 



where ^ 



For convenience, let us introduce the 
functions ft (i = 1,2,3,4) which relate the functions gt 
as the follows: 



91 

92 
93 



A, 

I 2 , 



-(mi - m 2 )\q\ 



-0-/2 + h 

mim 2 + u 1 Ld 2 + q z 



(cji+u 2 )|g1 
94 = ; h - .14 ■ 

From Eas.Hll|l). we may straightforward obtain the re- 
quirements fz = fi = and then from Eas. (|9ll()Jl and 
having the angular integration done, we obtain 

Mfi = m 1 +m 2 / 2 / -r^d\k\- — - — - 

it J \q] 2w 1 uj 2 {uji + io 2 ) 



+ 



|(mim 2 +ujilo 2 +q si )Qof 2 
(mi -m 2 ) 2 |g1|fc| j 

77117772 + E\E 2 +k 2 ' 



(15) 



M/2 = ^l±^t h 



mi -I- m2 
a s f \k\ 



d\k\—^ {(miW2 + m 2 wi)Qo/i 

2UJlUJ 2 I 

r jp \ i (16) 



I- j \q\ 

{Ei +E 2 )(w 1 +uj 2 )\q\\k\ 



here Ei = = y m ■ + P and the angular integra- 

tions of the equations have been carried out already, so 



iC+jfcj 

2M|fe| 



) (77 = 0, !,••■), the nth Legendre 
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functions of the second kind appear due to the precise 
BS kernel Eq.©. 

Eas. ((15l 116(1 are a complete set of coupled equations 
about /i and / 2 , and exactly equivalent to the instanta- 
neous BS equation. 

While from equations Eas. (|9HOI13J) i.e. the Breit equa- 
tion, with straightforward derivation, we 'directly' obtain 
the coupled equations: 



%d\k\ - m \ +m2 ; 



Mf a 
Mf 4 



(£i+£ 2 )(^i+^ 2 )|g1|fc| 
m 1 E 2 + m 2 Ei 

(ujx + u 2 )\q\Qifij , 

mi 

(mi - 



Qih 



-h 

■ m 2 
m 2 )f 3 ■ 



(18) 

(19) 
(20) 



M/i = {m 1 +m 2 )f 2 1 

7T 

• |(mim 2 + io\lo 2 + q 2 )Qof: 

(mi -m 2 ) 2 \q\\k\ 

+ — — ^-Wih 

m\m 2 + E\E 2 + k 

- (mi - m 2 )|q|<2i/3 j , 

M f2 = { ^±^h-mu 

mi + m 2 
ots_ f \k\ 

7T J \q\ "''"'2wiw 2 



Namely Eqs.|T7| d [H IHJ) are fully equivalent to 
Eci. (|12|) the Breit equation. It is easy to realize that 
if one sets fs = fi = 0, the equations Eas. lfTTl IThI H*9l 
I20|l for the Breit equation will 'return' to the equations 
Eas. l|15ll(j|l i.e. the instantaneous BS equation. 

The normalization condition for the solutions of the 
equations reads: 



(17) 



fd\q\ 



-{4g 2 [(m 2 - mi)\q\ 



.93 



d\k\— — |(miw 2 + to 2 wi)Qq/i 



(2tt) u>iuj 2 

+{oj 1 + Lu 2 )\q\gi + (miwg + m 2 wi)gi]} = 2M. (21) 



Table I. Expansion coefficients for the eigenfunctions (here means 0.000 • • •) 



{nl} 


WF 


Rio 


R20 


-R30 


i?40 


R50 


{10} 


ft 


0.707107 


7.59- 10~ 5 


4.34 ■ 10~ 5 


1.42 • 10~ 5 


4.26- 10~ 5 




h 


0.707107 


7.59- 10~ 5 


4.34 ■ 10~ 5 


1.42 • 10~ 5 


4.26- 10~ 5 




h 



















h 

















{20} 


h 


-7.98 • 10~ 5 


0.707107 


3.18- 10~ 4 


1.19- 10~ 4 


1.18 • 10" 4 




h 


-7.98 • 10- 5 


0.707107 


3.18- 10~ 4 


1.19- 10~ 4 


1.18 • 10" 4 




h 



















U 

















{30} 


h 


4.54 • 10~ 5 


3.19- 10~ 4 


-0.707106 


-5.2 • 10~ 4 


2.34- 10~ 4 




h 


4.54 • 10~ 5 


3.19- 10~ 4 


-0.707106 


-5.2 • 10~ 4 


2.34- 10~ 4 




h 



















h 

















{40} 


h 


-1.55 • 10~ 5 


-1.2 • 10~ 4 


-5.21 • 10~ 4 


0.707106 


7.66 • 10~ 4 




h 


-1.55 • 10~ 5 


-1.20- 10~ 4 


-5.21 • 10~ 4 


0.707106 


7.66 • 10~ 4 




h 



















h 

















{50} 


h 


-4.35 • 10~ 5 


-1.18- 10~ 4 


-2.33- 10~ 4 


-7.67- 10~ 4 


0.707106 




h 


-4.35 • 10~ 5 


-1.18- 10~ 4 


-2.33- 10~ 4 


-7.67- 10~ 4 


0.707106 




h 



















h 


















Now, let us solve the equations Eas. H17ll8ll9l20|) nu- 
merically by transforming the coupled equations into an 
eigenvalue problem of a matrix one, i.e., by expanding 
fi, i = 1, 2, 3,4 in terms of the bases of the exact S'-wave 
solutions of the Schrodinger equation (in momentum rep- 



resentation) R n i (|fc|)0: 
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Ak\ l c<l+l / n 2 [fc| 2 -l N 
(n 2 |fc| 2 + l)'+ 2 B_ '~ 1 Vn a |jf| a + l'' 

where CX^ix) is the Gegenbauer function, defined as the 
coefficient of h in the expansion of (1 — 2hx + h 2 )~ v 
in powers of h. To be practicable, as general cases, the 
expansion is truncated according to the request accuracy, 
so the present problem becomes an eigenvalue problem of 
a finite matrix. If we truncate the expansion up to j = 5 

if(k1)= £ 3?2r-*UI«1)> (23) 

j={nl} 

where j denotes the jth eigenvalue and cigcnfunction, n 



and I are the principal and angular quantum numbers: 
nl = 15, 25, 35, • • • (for _ states, only 5-wave states 
in the expansion are enough). By making the matrix for 
the four functions fi, (i — 1,2,3,4), corresponding to a 
(4 x 5) <g> (4 x 5) matrix to be diagonal, we finally obtain 
the results (eigenvalues and eigenfunctions accordingly). 
If ignoring the negative energy eigenvalue solutions which 
have similar physics meanings as those of Dirac equa- 
tions, we may organize the solutions into two types: type- 
A (the so-called 'trivial' solutions / 3 = / 4 = 0) and type- 
B (the so-called 'non-trivial' solutions fj, ^ 0, /4 ^ 0): 



Table II. Expansion coefficients for the eigenfunctions 



{nl} 


WF 


Rio 


-R20 


-R30 


R40 


-R50 


{10} 


h 


2.284- 10~ 3 


1.0679- 10- 3 


5.770 • 10- 4 


3.615 • 10~ 4 


2.407 • 10~ 4 


h 


2.262 • 10^3 


1.0577- 10- 3 


5.714- 10~ 4 


3.580 • 10- 4 


2.384- 10~ 4 


h 


-0.615535 


-0.284194 


-0.157829 


-0.101564 


-0.0715402 


h 


0.615535 


0.284194 


0.157829 


0.101564 


0.0715402 


{20} 


h 


4.169 • 10~ 4 


-3.866- 10~ 4 


-5.112- 10~ 4 


-3.862- 10~ 4 


-2.437- 10~ 4 


h 


4.129 • 10~ 4 


-3.829- 10~ 4 


-5.063- 10- 4 


-3.825- 10~ 4 


-2.414- 10~ 4 


h 


-0.320488 


0.306022 


0.4.3615 


0.309568 


0.211896 


h 


0.320488 


-0.306022 


-0.403615 


-0.309568 


-0.211896 


{30} 


h 


6.429 • 10~ 5 


-3.422- 10~ 4 


-2.655 • 10~ 5 


2.457- 10" 4 


2.267- 10~ 4 


h 


6.367- 10- 5 


-3.389- 10~ 4 


2.630 ■ 10- 5 


2.434 ■ 10~ 4 


2.245 • 10~ 4 


h 


-0.119399 


0.469465 


1.961 • 10~ 2 


-0.364445 


-0.363513 


U 


0.119399 


-0.469465 


-1.961 • 10~ 2 


0.364445 


0.363513 


{40} 


h 


-2.126- 10" 5 


1.1331 • 10~ 4 


-1.822- 10~ 4 


-3.558- 10" 5 


1.5455- 10" 4 


h 


-2.106- 10" 5 


1.1222 • 10~ 4 


-1.805- 10- 4 


-3.524- 10~ 5 


1.531 • 10~ 4 


h 


6.011 • 10~ 2 


-0.292407 


0.446233 


7.971 • 10~ 2 


-0.453219 


h 


-6.011 • 10" 2 


0.292407 


-0.446233 


-7.971 • 10~ 2 


0.453219 


{50} 


h 


1.786 • 10~ 7 


2.703- 10~ 5 


-5.647- 10~ 5 


9.173- 10~ 5 


-3.973- 10~ 5 


h 


1.769 • 10~ 7 


2.677- 10- 5 


-5.593 • 10~ 5 


9.085 ■ 10- 5 


-3.93438 


h 


2.300 • 10- 2 


-0.1403 


0.335672 


-0.504656 


0.335311 


h 


-2.300- 10- 2 


0.1403 


-0.335672 


0.504656 


-0.335311 



A. The so-called 'trivial' solutions 

The 'trivial'-eigenvalues (in unit eV) for (// ± e T ) sys- 
tems are 

-13.5410866, -3.3863112, -1.5048207 
-0.8467912, -0.5420610, 

The eigenfunctions corresponding to the eigenvalues 
have the expansion coefficients as those in Tabled. 

The over all renormalisation constants for the 'trivial' 
solutions are 45778.501, 45778.730, 45778.772, 45778.787 
45778.794, • ■ • accordingly. 

In fact, as expected and from Table I one may realize 
that the 'trivial solutions' of the Eas. (|17ll8ll9l20[> have 



fs = fi = exactly, and indeed they just correspond to 
the solutions of Eas. H15ll6|) . 

B. The so-called 'non-trivial' solutions 

The 'non-trivial'-eigenvalues (in unit eV) for (fj, e^) 
systems are 

-1.022015012- 10 6 , -1.021999529- 10 6 , 
-1.021998308- 10 6 , -1.021997962- 10 6 , 
-1.021997837- 10 6 , 

It means that the energy levels have so great binding 
energy almost as 2m e (two times of the electron mass). 
They are 'extra' to the instantaneous BS equation. It 
seems that in certain sense, the Breit equation is a 'bi- 
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product' of two Dirac equations, so there are extra solu- 
tions, which still 'fall' in the positive energy region for the 
binding systems, i.e. would correspond to those negative 
energy solutions if the lighter component were free. Thus 
the gap of the 'nontrivial solutions' from the 'trivial ones' 
can be so great as two times of the electron mass for the 
systems (/i ± e =F ). To understand the fact, we have done 
several exercises, and find that it is a general feature for 
the Breit equation that the 'trivial' and 'non-trivial' so- 
lutions have so deep a gap as the two times mass of the 
lighter component in the bound states. 

Since there is no evidence in experiments for the sys- 
tems (/i ± e =F ), so exactly to say, Breit equation is not 
physical. Only when one restrict oneself to consider the 
week binding solutions and ignore the (deep) 'nontrivial 
solutions', the Breit equation is meaningful. In literature, 
fortunately, people use Breit equation with a strong im- 
plication that only the very week binding solutions are 
pursued. Indeed there is no problem only when applying 
the Breit equation to the week binding spectrum stud- 
ies, whereas if applying it for 'complete' energy spectrum 
and/or complete set of the wave functions to consider 
relativistic correction effects etc, it will cause problem. 
However the instantaneous BS equation does not have 
the problem at all. 

The eigenfunctions corresponding to the eigenvalues 
have the expansion coefficients as those in Table. II. 

The over all renormalisation constants for the 'non- 
trivial' solutions are 

1.240240535 ■ 10 7 , 3.652200818 • 10 7 , 6.728317282 • 10 7 , 
1.215161253 • 10 8 , 2.757689877- 10 8 , • ■ • accordingly. 

To test the 'stability' of the numerical solutions, we 
also try to truncate the expansion up to j — 10, and 
obtain the eigenvalues for 'non-trivial' and 'trivial' as the 
follows: 

-1.022015269 • 10 6 , -1.021999737 • 10 6 , 
-1.021998500 • 10 6 , -1.021998133 • 10 6 , 
-1.021997979 • 10 6 , -1.021997900 • 10 6 , 



-1.021997856 • 10 6 , -1.021997830 • 10 6 , 
-1.021997815 • 10 6 , -1.021997807 • 10 6 , 



-13.5410866, -3.3863112, -1.5048207, 
-0.8467912, -0.5420610, -0.3764843, 
-0.2764580, -0.2116343, -0.1672546, 
-0.1353994, 

From the eigenvalues obtained by a different truncation, 
one may see the low-laying numerical solutions for j = 5 
are quite accurate for our observation. The table for 
j = 10 cigenfunction coefficients is too big to present 
here, thus we omit it. 

From the example, we have learnt that the precise 
numerical solutions of the equations Eas. H17ll8ll9l20|) 
have confirmed the arguments on the Eas. (|13fl about the 
homogeneous equations. Generally, the Breit equation 
Eq. (|12|) contains not only all the solutions of the instan- 
taneous BS equation, but also the un-physical solutions, 
the 'non-trivial' solutions, although we have the lessen 
only for J p — 0~ states of the (/i ± e T ) systems. There- 
fore, we can conclude that the Breit equation Ea. l|12|) 
is not equivalent to the instantaneous BS equation, and 
the Breit equation Eci. (|12fl should be abandoned for all 
kinds of instantaneous interacting binding systems, if one 
wishes to apply it to finding exact and complete solutions 
and/or uses it for computing the relativistic correction ef- 
fects without restriction. In contrary, the instantaneous 
BS equation does not have such un-physical solutions 
for the instantaneous interacting binding systems at all, 
therefore one may take liberties with applying the BS 
equation to various problems of the systems. 
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